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ABSTRACT 


I  he  results  of  a  study  of  vs  m^  characteristics 
ot  earthquakes  in  the  liastern  Himalayan  region  arc  given 
in  this  report.  It  is  shown  that  in  this  region  some 
earthquakes  occur  which  have  vs  m^  characteristics 
similar  to  explosions,  exhibiting  low  surface  wave 
magnitudes  relative  to  body  wave  magnitudes,  vhen  seen 
at  the  reporting  stations  available  to  this  study.  The 
application  of  station  corrections  does  not  change  the 
general  distribution  and  spread  of  points  in  the  M. 
vs  m^  plane,  and  therefore  it  is  unlikely  that  station 
(or  path)  effects  are  the  source  of  the  anomalies,  bocal 
depths  of  most  of  the  events  studied  are  shallow  or 
normal  and  can  also  be  ruled  out  as  causes  of  low 
surface  wave  magnitudes.  The  geographical  distribution 
of  anomalous  events  correlates  with  various  prominent 
geological  features  and  probably  reflects  the  distri¬ 
bution  of  tectonic  stress  in  the  area  studied.  The 
existence  of  th^se  anomalous  events  in  certain  areas 
of  the  world  can  seriously  decrease  the  effectiveness 
of  the  Mg  vs  criterion  in  discriminating  between 
earthquakes  and  explosions.  Detection  of  the  Rayleigh 
wave  from  these  events  is  so  difficult  that  further 
understanding  of  their  mechanism  is  dependent  on  an 
improved  monitoring  capability  with  good  azimuthal 
coverage  „ 
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•  NTROIMK’T  ION 


This  report  describes  .1  study  of  Jl  vs  a.  tb.1r.1c- 

js»  n 

teristics  of  ea  tliquakes  in  a  United  region  comprising 
tht  eastern  Itiiaa layas ,  and  parts  of  Assam  and  Tibet. 

Ibis  study  i*i  preliminary  and  is  intended  to  be  con¬ 
tinued  and  expanded  in  scope. 

v  s  charat  terist  ics  are  general  ly  thought  to  be  out 
of  the  most  effective  seismic  discr lain.mts  between  earth¬ 
quakes  and  explosions.  I.arthquakes ,  in  general,  excite 
I  tug  period  surface  waves  much  aore  efficiently  than 
do  nut  1  ear  explosions.  The  relative  excitation  of  short  - 
period  body  waves  and  long-period  surface  waves  by 
earthquakes  is  likely  to  be  determined  by  the  charac¬ 
teristics  of  the  source,  especially  the  stress  in  the 
source  region  ( It  rune  IDbH;  hvss  and  Brunt*  1908;  Kyss 
I9?0a,  1 S* 7 III* ;  lh.it  her  1972;  lhatcher  and  Itrune  I «»7 1 ; 

Hanks  11*71.)  The  location  of  earthquake  regions  amt 
their  source  characteristics  are  no  thought  to  he 
determined  hy  the  distribution  ol  lithospheric  plate 
hoiindaries  and  their  relative  notion  (1  sacks,  Oliver, 
and  Sykes  ll>(»8;  l.el'iclion  ll>68). 

It  can  he  expected,  therefore,  that  the  source 
character  is t *cs ,  in  particular,  the  stress  present  in 
the  source  region  will  vary  cons iderab I y ,  depending  on 
the  particular  region  considered.  Decent  research  has 
shown  that  there  Is  a  wide  variation  in  the  relative 
amplitude  of  short-period  body  waves  and  long-period 
surface  waves  from  earthquakes  which  cannot  he  explained 
by  differences  in  propagation  paths  or  radiation  patterns 
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(Ihatcher  1 1> 7 2 ;  Ihatcher  ;uul  Bruno  1971;  Marshall  and 
Basham  1972).  Such  variations  arc  of  considerable 
interest  in  the  problem  ol  detection  and  identification 
ol  nuclear  explosions,  since  high -stress  earthquakes 
could  be  mistaken  for  nuclear  explosions.  On  the  other 
hand,  nuclear  explosions  could  be  effectively  hidden  in 
hijih-strcss  regions.  I  rom  a  long  range  point  of  view  it 
is  also  very  important  to  map  the  distribution  of  stress, 
since  it  contains  information  about  do  format iona  1  pro¬ 
cesses  currently  in  progress  within  the  earth. 
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TECHNICAL  DISCUSSION 


Data 


Ihc  region  studied  is  bounded  by  latitudes  27°N 
and  34°N  and  longitudes  92°I;  and  1 0 0 0 1; ,  and  is  shown 
in  ligurc  1  along  with  the  stations  used.  Geographically 
tills  region  covers  the  eastern-most  part  of  the  Himalaya 
range  and  part  of  the  north-south  trending  mountain 
ranges  which  continue  southward  to  Burma.  Toward  the 
soutii,  part  of  Assam  is  included,  and  in  the  north  a 
portion  of  the  Tibetan  Ulatcau.  The  available  geological 
information  is  not  sufficient  to  describe  the  region 
accurately,  as  would  be  desirable.  (Gansscr  1964; 

Smirnov  1964.)  A  tectonic  map  of  the  region  has  been 
published  by  the  Geologic  Institute  of  the  Academy  cf 
Sciences  of  the  USSR.  The  Tibetan  portion  of  the  map 
is  based  on  recent  Chinese  works  in  the  area.  The  map 
is  reproduced  in  Figure  2. 

Tapes  containing  NOS  cpicentral  data  have  been 
searched  for  earthquakes  occurring  within  this  region 
between  1960  and  1971,  Approximately  one-hundred  events 
were  found.  Data  for  67  events  was  ordered.  Most  of  the 
remaining  events  were  not  ordered  either  because  NOS 
did  not  compute  a  body  wave  magnitude,  or  because  their 
magnitude  was  less  than  4.0.  The  67  events  are  listed 
in  Table  IA„  Coordinates  and  elevations  for  the  sta¬ 
tions  used  arc  give/,  in  Table  IB.  The  same  table  also 
shows  the  average  distances  and  azimuths  from  the  center 
of  the  region  studied  to  all  the  stations.  Eleven  of 
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the  67  events  were  not  used  in  the  analysis,  however, 
due  to  lack  of  or  poor  quality  of  data.  Two  additional 
events  were  not  analyzed  because  they  were  mixed  with 
other  events.  The  total  13  events  are  indicated  by 
daggers  in  Table  IA.  Most  of  the  events  have  shallow 
or  normal  focal  depths,  and  are  scattered  over  all 
parts  of  the  region  but  there  are  the  following  main 
areas  of  seismicity:  the  Himalayan  front,  the  eastern 
end  of  the  plain  of  Assam  and  mountain  ranges  trending 
southeast  towards  Burma,  and  a  region  at  the  junction 
of  the  mountain  ranges  at  around  30°N  and  95°  1*  which 
exhibits  swarmlike  earthquake  activity.  There  is 
also  a  concentration  of  events  around  33°N  and  93°li. 
figure  lb  from  Yanshin  (1966)  shows  that  30°N  9 5 ° Ii  to  be 
close  to  a  major  fault  which  bynr  Sykes  (personal 
communication)  suggests  may  be  the  India-Asia  continental 
block  suture.  The  rest  of  the  activity  is  scattered; 
there  are  some  possible  zones  of  activity  and  inactivity 
whic.i  cannot  be  seen  clearly  because  of  the  paucity  of 
events.  As  pointed  out  by  Fitch  (1970),  the  region  of 
Tibet  and  Western  China  is  Very  complex  and  the  details 
of  regional  plate  motions  have  not  been  worked  out  yet. 
The  situation  is  different  for  the  Himalayas  and  the 
crustal  block  of  India.  The  recognition  that  India  is 
thrust  under  the  Tibetan  plateau  predates  the  notion 
of  plate  tectonics  by  many  years.  The  majority  of  events 
used  in  this  study  were  listed  by  NOS  as  shallow  or 
normal  depth  earthquakes.  These  depths  were  presumably 
determined  from  travel  times  and  pP  phases.  Only  one 
earthquake  has  a  depth  exceeding  100  km,  and  only  a 
few  are  deeper  than  50  km.  Depth,  therefore  as  we  shall 
see  in  more  detail  later,  cannot  be  an  explanation  for 
the  weak  surface  waves  from  some  of  the  events. 
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For  all  the  events,  film  chips  were  obtained  from 
the  National  Ocean  Survey  for  the  following  stations: 
Shillong  (Sill.),  New  Delhi  ( N I )  I )  ,  Kabul  (KBL)  ,  Quetta 
(QUIi)  ,  Lahore  (LAII)  ,  Poona  (.’00),  and  OhiangMai  (0110). 
These  stations  are  equipped  with  the  standard  WWSSN 
instrumentation.  The  period  of  the  Id’  seismograph 
pendulum  was  changed  during  the  year  1965  from  30  seconds 
to  15  seconds.  This  improved  the  detection  ol  short 
period  (10-15  second)  surface  waves  which  lor  some  small 
events  are  the  most  prominent.  Since  most  events  of 
interest  were  small  in  magnitude,  it  was  decided  not  to 

order  data  for  more  distant  stations. 

c 

Unfortunately  there  is  not  enough  data  available 
to  establish  regional  body  wave  magnitude  formula  in 
the  areas  studied,  as  was  done  by  Lvernden  (1967)  for 
the  United  States.  The  complexity  of  this  region  makes 
it  likely  that  several  different  crustal  models  would 
be  necessary  to  explain  crustal  phases  at  all  stations. 

The  coverage  of  the  stations  is  not  sufficient  to 
establish  such  models.  Attempts  to  find  systematic 
patterns  of  l’  phases  at  various  stations  and  thus 
separate  them  from  depth  phases  yielded  no  conclusive 
results.  Search  for  pi’  phases  also  did  not  give  consistent 
results . 

Computations 

Magnitudes  were  computed  for  all  events  in  the 
sample.  Short-period  body  wave  magnitudes  (mb)  were 
computed  with  the  conventional  Gutenberg  formula,  in 
the  absence  of  better  body  wave  magnitude  formulas  for 
the  region.  This  undoubtedly  contributes  to  the  scatter 
of  mb  values. 
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To  evaluate  the  network  bias  in  we  selected  the 
six  ouc  oi  the  original  67  events  Tor  which  the  NOS 
Earthquake  Data  Reports  show  eight  or  more  stations  used 
in  the  calculation  for  m^.  For  five  of  the  six  events 
we  also  were  able  to  calculate  an  value,  (Station 
corrections  were  calculated  and  used  to  remove  any 
station  bias  with  respect  to  the  network,  this  would  not 
of  course  remove  any  network  bias,)  On  the  average  the 
NOS  magnitudes  were  larger  by  0.2  ±  0,1  (one  standard 
deviation  of  the  mean).  This  value  becomes  0,1  ±  0.15  if 
stations  for  which  A>90°  or  A<20°  are  excluded  from  the 
NOS  magnitude  calculations  but  all  other  stations 
reporting  magnitudes  are  included.  If  our  magnitudes 
were  corrected  to  remove  this  estimated  bias,  the  earth¬ 
quakes  would  look  more  explosion-like. 


Had  we  used  Veith  and  Clawson's  (1972)  mb  curve 
instead  of  Gutenberg's  curve,  our  network  bias  would 
have  been  reduced  by  0.05  m^ . 

The  surface  wave  magnitudes  were  computed  by  three 
methods,  (1)  using  the  (modified)  Prague  formula: 


Ms  =  log  (A/T)  +  1.66  logA  -  0.18 

where  A  is  the  maximum  peak-to-peak  amplitude  in  milli¬ 
microns,  (2)  the  method  introduced  by  Marshall  and  Basham 
(1972)  for  near  distances: 

Ms  =  log  A  +  B(T)  +  C (A) 

where  A  is  the  maximum  zero-to-peak  amplitude  in  milli¬ 
microns,  with  a  depth  correction  given  by  the  equation 


tfs  -  Ms  ♦  .008h(km) 


(U 


(2) 


(2a) 
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where  h  is  determined  from  values  taken  at  different 
peiiods,  and  (3)  a  formula  proposed  by  von  Scggcrn  (1970) 
for  near  distances: 

M  =  leg  (A/T)  +  1.16  1 o  g  A  +  0.7  4 


where  A  is  the  maximum  pcak-to-peak  amplitude  of  the 
surface  wave  train  (in  millimicrons),  T  is  the  wave 
period,  and  A  is  the  epicentral  distance  in  degrees. 
B(T)  and  ('(A)  are  the  period  and  distance  correction 
factors  which  are  tabulated  by  Marshall  and  Basham.  In 
our  calculations  we  used  their  tables  computed  for 
Central  Asian  surface  wave  dispersion  characteristics. 
Other  equivalent  approaches  to  the  regionalization  of 
magnitude  calculations  with  respect  to  dispersion 
characteristics  were  reported  by  Alewine  (1972)  and 
Basham  (1971).  Of  the  methods  used, only  Marshal]  and 
Basham's  method  attempts  to  correct  for  deptn. 

Rcsul ts 


The  plots  of  Mg  versus  in^  are  given  in  figures  3, 

4  and  3.  Those  magnitudes  were  computed  by  averaging  the 
magnitudes  at  all  available  stations,  figure  3  shows  the 
surface  wave  magnitude  computed  from  the  Prague  formula. 
Since  the  stations  used  arc  at  small  epicentral  distances 
for  the  area  studied,  this  formula  tends  to  give  low  Ms 

estimates.  Of  course  if  explosions  were  recorded  at  the 

same  distance,  thev  also  would  have  low  M  values  and 

s 

the  separation  would  be  unaffected,  figures  4  and  5 
respectively  show  the  magnitudes  computed  by  the  methods 


(3) 
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given  by  Marshall  am!  Basham  (  1 (.  V .! )  with  depth 
corrections,  and  von  Seggern  ( 1  J 7 0 ) .  ihe.se  formulas 
were  designed  to  eliminate  the  depend  mice  o  1'  surface 
wave  magnitudes  on  epic ‘Mitral  <.i  i  nance  at  short 
distances  ( A  2  S ' '  J  and  are  therefore  more  ippJicable 
the  problem  investigated.  1  rend  lines  obtained  by 
Marshall  and  Basham  (  1 U 7 2 )  using  data  reported  by  Capon 
et  al.  (10(>7)  for  M  vs  m^  dependence  of  earthquakes  and 
explosions  in  Central  \sia  are  superimposed  on  all  these 
figures,  the  magnitudes  given  are  scattered  between  the 
two  lines,  aid  although  most  of  them  are  above  the  explo¬ 
sion  lines  some  of  them  are  fairly  close  to  it.  Magnitudes 
computed  by  Marshall  and  Basham's  method  are  the  least 
scattered  of  all  three  : ets  of  magnitudes  shewn. 

For  til  events  for  which  body  waves  could  be  detected 
at  any  station,  surface  waves  could  also  be  detected  at 
SIIL  when  Sllh  was  operational.  In  general,  for  those  events 
for  which  surface  waves  could  he  detected  only  at  SIIL, 
the  noise  was  normal  at  the  other  stations,  suggesting 
that  the  signals  truly  were  small  and  that  the  small 
amplitude  at  Sllh  is  not  the  effect  only  of  a  radiation 
pattern. 

Table  IV  shows  all  the  computed  magnitudes  together 
with  averages. 

In  order  to  depict  the  geographical  distribution 
of  the  magnitude  characteristics,  the  M  vs  m^  plane  in 
all  cases  was  subdivided  into  four  parts  by  the  lines 
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=  mb  ' 
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=  ,ub  " 
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(4) 

(5) 

(6) 


In  order  to  avoid  cluttering  the  figures,  these 
lines  are  not  shown  on  the  M  vs  m(  plots,  but  events 
1 rom  the  four  parts  of  the  vs  planes  are  plotted 
with  different  symbols  on  a  map  having  equal  latitude 
and  longitude  intervals,  fince  the  region  studied  is 
relatively  close  to  the  equator,  this  does  not  cause 
a  great  amount  of  areal  distortion.  The  squares  on 
the  plot  indicate  events  which  are  the  most  expios ionlike , 
and  fall  below  the  line  defined  by  equation  (6),  diamond 
symbols  denote  events  in  the  next  part  of  the  M  vs  m, 
plane  above  this,  crosses  denote  events  falling  into  the 
next  division  above,  and  finally,  events  which  ar'>  above 
the  line  defined  in  equation  (4)  are  denoted  by  vertical 
straight  lines.  figures  24a  through  24e  show  some  examples 
of  seismograms  which  show  the  obvious  differences  in 
the  Rayleigh  (and  Love)  wave  excitation  of  events  with 
similar  body  wave  magnitudes.  The  event  on  the  left  in 
figures  24a, b,c  show  large  Rayleigh  waves  and  a  large 
hove  wave  pulse  which  is  absent  in  the  event  on  the 
right.  S  and  hove  wave  excitation  is  another  matter  which 
should  be  investigated  further.  However  we  may  say  here 
that  for  none  of  the  events  represented  by  a  square  in 
figures  (> ,  7,  and  8  could  S  waves  or  Love  waves  greater 
than  the  Rayleigh  wave  be  seen. 

figures  (>,  7,  and  8  show  the  plots  obtained  by 
using  various  magnitude  formulas.  Although  the  classi¬ 
fication  of  some  events  changes  depending  on  the 
magnitude  formulas  used,  all  plots  show  essentially  the 
same  geographical  distribution  pattern,  for  this  type 
of  figure  a  transparent  overlay  is  included  in  the  back 
pocket  to  facilitate  geographical  orientation. 


The  most  striking  feature  of  the  map  is  the  high 
concentration  of  explos Lonlikc  events  in  an  elongated 
reg  i  on  centered  at  about  30°N  and  95°!:.  The  earthquakes 
shown  in  this  region  occurred  mainly  in  two  sequences 
of  events  with  body  wave  magnitudes  ranging  from  4.5  to  5, 
one  of  which  occurred  in  .June  to  August  1968,  the  other 
in  .June  to  September  1969. 

Besides  this  region  there  are  two  clear  regions 
which  are  also  characterized  by  low  M  values  relative 
to  ii»h  values.  One  is  the  frontal  region  of  the  Himalayas 
also  described  by  Marshall  and  Basham  (1972),  which  shows 
such  events  intermingled  with  normal  earthquakes.  Another 
event  was  reported  recently  in  this  region,  on  October  24, 
1971  at  28. 2N  and  87. 21;  with  m})  (NOS)  =  5.1  and  M.  2  5 
at  LASA,  NORSAR  and  (111(1 .  The  other  occupies  the  eastern 
end  ol  Assam  and  north-south  trending  mountain  ranges 
which  join  the  Himalayas  an  1  extend  south  towards  Burma. 

I  he  anomalous  region  at  30oN  and  9 5 ° I :  mentioned  above  is 
situated  at  -lie  eastern  end  of  the  high  mountains  belonging 
to  the  Himalayas,  where  the  trend  of  the  mountains  turns 
towards  the  south.  And  a  deep  fault  terminates  against 
a  (Irani  to  id  body.  (See  Figure  2.) 

The  observed  magnitude  characteristics  are  thus 
seen  to  be  correlated  with  geological  features,  probably 
re!  looting  distribution  of  tectonic  stress.  Other  trends 
may  also  be  present  but  they  arc  not  consistent  on  the 
various  plots  presented  here,  and  the  number  of  cvcn-s 
is  not  sufficient  to  obtain  a  clear  pattern.  The  interior 
region  of  Tibet  seems  to  be  dominated  by  shallow  earth¬ 
quakes  displaying  normal  earthquake  M  vs  m^  characteristics. 
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The  Teast  squares  lines  fitted  to  M  vs  in,  values 

s  L) 

of  earthquakes  and  explosions  in  different  regions  of  the 
earth  have  had  slopes  varying  between  1.0  and  2.0  with 
most  of  the  reported  values  being  slightly  above  1.  Since 
in  the  previous  figures  the  vs  plane  was  subdivided 
with  lines  having  a  slope  o  1.0,  it  is  necessary  to  test 


whether  lines  with  a  different  slope  would  change  the 
general  picture,  figures  9,  10,  and  11  show  results 
obtained  when  the  vs  m,  plane  was  subdivided  with  the 

1  incs 

M 

s 

=  1,5  m,  -  5.6 

(7) 

M 

s 

=  1.5  mb  -  3.1 

fHJ 

M 

s 

=  1.5  mb  -  2.6. 

(9) 

l'hc  symbols  were  designated  as  before  relative 
to  these  new  lines.  The  figures  show  that  although  the 
classification  ol  many  events  changed,  the  general 
pattern  is  the  same.  That  is,  explosionlike  events 
falling  into  the  lowermost  regions  in  the  M  vs  mb  plane 
are  still  concentrated  in  the  same  geographical  regions. 
I  here  1  ore ,  it  is  concluded  that  the  slope  of  lines  for 
the  given  data  set  does  not  affect  the  geographical 
distribution  of  various  types  of  events. 

The  average  magnitude  calculations  presented  above 
arc  influenced  by  path  and  station  effects.  The  station 
ellects  cause  bias  due  to  incompleteness  of  data,  since 
a  varying  number  of  stations  arc  available  for  indi¬ 
vidual  events.  In  order  to  estimate  the  station  effects 
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the  difference  between  magnitudes  at  various  stations 

were  averaged.  Table  II  shows  mean  magnitude  differences 

using  Quetta  (QUIi)  as  a  reference  station.  The  standard 

deviation  o'  the  differences  from  the  mean  are  also 

shown.  Most  differences  were  estimated  using  about 

thirty  events,  except  those  involving  KBL  which  involved 

only  about  (i-7  events.  The  table  shows  that  SI  IL ,  CMC 

and  KBL  tend  to  give  considerably  higher  nij  estimates 

than  the  rest  of  the  stations,  ,.ml  the  same  stations 

tend  to  give  low  M  values.  Therefore  earthquakes  whose 

M  vs  m^  values  arc  based  on  only  these  stations  will 

appear  more  cxplos ion  1  ike .  M.  corrections  for  Marshall 

and  Basham's  M  formula  arc  the  smallest  of  the  three 
s 

formulas  considered,  indicating  that  this  one  is  the 
most  appropriate  for  determining  M  for  this  region. 

The  Prague  formula  gives  too  small  M.  values  at  short 
distances,  while  von  Seggern's  formula,  based  on  NTS 
and  Western  U.S.,  tends  to  overcorrect  and  gives  too 
high  M  values  at  short  distances. 

It  was  decided  to  define  station  corrections  in 
such  a  way  that  the  mean  of  the  corrections  is  zero, 
that  is,  if  all  stations  in  the  data  group  arc  avail¬ 
able,  the  mean  magnitude  remains  unchanged.  This  method 
would  make  corrections  for  data  with  incomplete  sets  of 
stations  hut  would  not  change  the  bias  of  the  network 
used.  The  corrections  thus  defined  for  M  turn  out  to 
he  fairly  small,  hut  some  of  the  corrections  for  m^  are 
large.  Nevertheless,  the  bias  introduced  by  the  small 
network  is  judged  to  he  small,  since  over  all  average  in^ 
values  are  very  close  to  the  NOS  values,  which  arc  typi¬ 
cally  based  on  observations  at  3-5  telcscismie  stations. 
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liven  if  the  station  magnitude  differences  between  KBL 
and  the  rest  of  the  stations  were  disregarded,  which 
could  be  just.fied  on  the  basis  of  too  few  observations 
for  this  station,  the  average  values  would  not  decrease 
by  more  than  0.2  magnitude  units. 

lable  V  shows  the  corrected  magnitudes.  Marshall 
and  Basham's  depth  correctiqn  lias  been  applied  by  using 
the  NOS  reported  deptli  (33  km  for  normal  events)  in 
formula  2a.  Measurement  at  the  different  Rayleigh  wave 
periods  required  to  estimate  h  was  often  impossible 
because  of  the  pulse-like  character  of  the  close-in 
Rayleigh  waves  which  were  the  only  ones  seen  for  the 
anomalous  events.  The  corrections  reduce  the  average 
standard  deviation  of  m^  at  individual  stations  around 
the  event  mean  values  to  one-half  of  the  original  value 
and  are  thus  fairly  effective,  figures  12,  13,  and.  16 
are  the  corresponding  Mg vs  plots. 

Since  the  Marshall  and  Basham  method  uses  a  depth 
correction  the  average  value  determined  by  this 
method  is  relatively  high.  However,  Marshall  and 
Basham  (1972)  found  that  application  of  the  depth 
correction  does  not  cause  a  shift  in  the  explosions  but 
only  in  the  earthquakes,  thereby  increasing  the 
separation  between  the  two  populations. 

Further  light  is  thrown  on  the  effect  of  depth  on 
vs  m^  by  Figure  14  where  depth  corrections  have  not  been 
applied,  but  where  the  events  of  greater  than  60  km  depth 
have  been  indicated  by  squares;  and  the  events  of  unknown 
depth  by  filled  upright  triangles.  We  see  that  there 
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arc  many  shallow  events  near  the  explosion  population 
mean.  Of  course  if  by  using  digital  recording  the 
spectral  depth  estimate  could  he  made,  and  if,  as 
Marshall  and  Basham  assert,  the  explosion  population 
docs  not  move  and  if  the  anomalous  events  do  move  up 
on  the  Ms  vs  m^  plot,  then  Figure  14  would  he  much  like 
Figure  13  and  show  better  separation.  In  Figure  15  the 
Marshall  and  Basham  values  are  plotted  against  NOS 
in.  values.  This  figure  may  be  compared  with  Figure  13. 

In  Figure  13  there  arc  no  events  for  which  M  <  -  1,5, 

although  three  arc  very  close.  Thus  use  of  NOS  m^  values 
results  in  events  being  more  anomalous  than  does  use  of 
the  SDL  m^  values.  This  is  in  accord  with  the  bias 
calculations  in  the  Computations  section.  As  an  example 
of  the  dangers  of  using  NOS  magnitudes,  consider  the 
point  in  the  stippled  region  with  in^  =  5.5.  This  in^ 
value  comes  from  one  observation  at  SIIL,  a  distance 
of  7.1°.  The  SDL  magnitude,  with  station  corrections, 
is  4„32.  This  problem  is  also  present  in  the  work  of 
Landers  (1972)  where  his  NOS  magnitudes  for  Tibetan 
events  average  0.25  m^  higher  than  ours.  Because  he 
did  not  apply  depth  corrections,  his  Mg  values  are  on 
the  average  0.23  Mg  lower.  Thus  his  results  arc  about 
0.5  magnitude  units  "more  anomalous"  than  ours. 

Figures  17,  18,  19,  and  20  show  the  geographical 
distribution  of  various  types  of  events  as  defined  by 
lines  of  unit  slope  described  above  after  corrections 
have  been  applied.  Although  classification  of  individual 
events  changes  somewhat,  the  general  pattern  did  not. 
Although  in  Figure  18  no  events  with  observed  Rayleigh 
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waves  remain  ior  which  Ms  <  ra,,  -  1.5,  three  of  them  arc- 
very  close  to  the  line. 

Station  corrections  were  also  computed  hy  the  joint 
magnitude  determination  (JMH)  method  (von  Seggcrn,  1972) 
using  only  events  recorded  by  at  least  three  stations. 
The  corrections  and  the  slopes  and  intercepts  ol  the 
least  squares  linear  fit  to  points  in  the  log  A  -  log  \ 
plane  are  given  in  Table  Ill  (A  is  the  amplitude  of  1* 
and  Rayleigh  waves  respectively.)  Figures  21a  and  21b 
show  plots  of  log  A  before  and  after  the  corrections  for 
1’  were  applied.  A  visible  decrease  of  scatter  results  il 
station  corrections  are  used. 


Figures  22a  and  22b  show  a  similar  plot  for  Rayleigh 
wave  amplitudes.  In  this  case  the  decrease  of  scatter  due 
to  the  application  of  corrections  is  negligible. 

Figure  23  shews  a  plot  of  event  factors  (!:s  for  surface 
waves  versus  F,  for  I’  waves)  which  is  analogous  to  the 
M  -  plots.  Triangles  denote  events  whose  surface 
waves  were  recorded  at  ’css  than  three  stations,  and 


these  events  were  excluded  for  the  .JMI)  calculations. 


Since  many  of  these  events  have  anomalously  low  surface 
wave  amplitude?,  and  arc  therefore  of  prime  interest 
in  this  study,  their  relative  event  factors  were 


computed  subsequently  using  the  exponential  decay  factor 
and  station  corrections  derived  for  the  more  complete 
data  set.  It  was  found  that  the  relative  positions  and 
the  scatter  of  events  in  the  \\.  -  l:b  plane  is  similar 
to  those  in  the  Hgvs  i«b  plots,  and  the  geographical 
distribution  of  explosion  and  earthquake  type  events  is 


unul tc  red . 
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It  is  interesting  to  speculate  about  the  reasons 
lor  the  explos i on  1  ike  character  of  many  events  of  this 
region.  Besides  the  high  stress  at  the  source  it  is 
possible  that  in  the  surrounding  regions,  especially  in 
the  thick,  high-Q  crust  of  Tibet,  the  propagation  of 
high-frequency  I*  waves  is  especially  effective,  and 
this  could  cause  the  vs  values  to  be  more  explosion' 
like.  This  supposition  is  contradicted,  however  by  the 
dominant  periods  oi  the  1’  waves,  which  range  between 
0.5  -  2  seconds  with  an  average  around  1.2  seconds, 
while  P  waves  in  regions  of  high  Q  show  much  higher 
frequencies.  (Isacks  et  al.  1968.)  That  the  anomalous 
events  are  indeed  earthquakes  and  not  explosions  is 
indicated  by  the  fact  that  many  arrivals  on  the  short 
period  instruments  are  dilatations.  A  study  of  these 
ana  other  phases  is  continuing.  It  will  bo  difficult 
however  to  do  a  truly  detinnivo  study  until  i„,provcd 
stations  and  arrays  arc  available  in  the  area. 

It  must  also  bo  noted  that  the  maximum  Rayleigh  wave 
amplitudes  of  many  small  events  occurred  at  very  short 
periods  (8  -  IS  seconds),  for  deep  events  the  dominant 
period  should  have  been  larger.  The  dependence  of  surface 
wave  excitation  on  depth  is  more  critical  at  short 
periods  and  bettor  depth  determination  could  improve  M 
values.  Whether  this  would  change  the  character  of  the" 
anomalous  events  is  not  clear  at  this  point. 

Wc  should  emphasize  that  the  depth  corrections  used 
in  this  report  were  made  using  NOS  depths,  not  M  measure¬ 
ments  as  done  by  Marshall  and  Sasham.  Thus  discrimination 
cou.d  only  he  aided  by  this  technique.  Wc  were  unable  to 
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measure  a  sufficient  number  of  periods  for  any  but  a  few 
events  to  determine  a  depth  independent  of  NOS.  Strictly 
following  Marshall  and  Basham's  procedure  would  have 
resulted  therefore  in  almost  no  depth  corrections  being 

applied. 
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CONCLUSIONS 


Some  earthquakes  in  the  cistern  Himalaya  region 
exhibit  anomalously  low  Mg  values  relative  to  values. 
The  geographical  distribution  of  such  events  shows  a 
pattern  which  correlates  with  various  geological 
features,  the  frontal  region  of  Himalayan  mountain 
ranges  and  the  sudden  change  of  the  trend  of  mountain 
ranges  at  30°N  and  95CE.  Application  of  station  or 
depth  corrections  does  not  change  the  general  geo¬ 
graphical  pattern  of  anomalous  events.  We  tentatively 
conclude  that  the  anomalous  Mg  vs  m^  character  is  due 
to  high  stresses  in  the  source  region.  The  absolute 
values  of  Mg  and  m^  are  such  that  they  resemble  explo¬ 
sions  detonated  in  other  regions  of  the  earth. 
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EARTHQUAKE  AND  EXPLOSION 
MEAN  TREND  LINES  FROM 
MARSHALL  AND  BASHAM  (1972) 
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text.  Prague  formula  was  used  for  Ms . 


NORTH  LATITUDE 


37.0 


35.5 

1  1  1 

_  n  Ms<mjj-1.5 

o  m|j-1.5<Ms<mjj-1 

I  1  1 

x  m|,-1<IVIs<m|j-.5 
'  mb-.5<Ms 

1 

34.0 

1 

1 

— 

32.5 

,  * 

X 

— 

31.0 

1  l 

_ _  x  l 

11-0 

6-  x 

— 

29.5 

1 

—  0 

X 

II 

1  X 

1 

28.0 

o 

s 

* 

J _ 

1 

1 

O 

26.5 

— 

1 

25.0 

__J _ 1 _ 1 _ 

1  1  1 

J _ 

90.0  91.5  93.0  94.5  96.0  97.5  99.0  100.5  102.0 

EAST  LONGITUDE 


Figure  7.  Geographical  distribution  of  various  Ms  vs  types. 
Definition  of  types  is  determined  by  equations  4,  5,  and  6  in 
text.  Marshall  and  Basham  method  was  used  lor  Mh . 
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I  iguro  K.  (Jeograph Leal  distribution  of  various  Ms  vs 
Definition  of  types  is  determined  by  equations  4,  5, 
text.  Von  Seggern's  formula  was  used  for  . 
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I  igure  9,  Geographical  distribution  of  various  Ms  vs  mp  types 
Definition  of  type  is  determined  by  equations  7,  8,  and  9  in 
text.  Prague  formula  was  used  for  Ms. 
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I'ij’ure  10.  (leographical  distribution  of  various  Ms  vs  types 
Definition  of  type  is  determined  by  equations  7,  8,  and  l.)  in 
text.  Marshall  and  Basham  method  used  for  Me. 
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I  i  iiiir-o  11.  (Icograplii  cnl  distribution  of  various  Ms  vs  inn 
Del  i  ii  it  ion  ol  type  is  determined  bv  equations  7,  8  and 
text.  Von  Severn’s  formulas  was  used  for  M  . 
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iigUit.  |5.  Mm  computed  with  the  method  oi'  Marshall  and  Ha  sham 
I  Vi  i  t  h  depth  corrections)  plotted  with  NOS  mb  values.  Stippled 
region  ‘'hows  area  occupied  by  explosions  in  the  paper  by 
Marshall  and  Hasham  ( 1 D  7  2 ) .  Bars  give  °nc  standard  deviation 
of  the  mean  for  selected  events.  (♦  <?/.!) 
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Figure  18.  Geographical  distribution  of  various  Ms  vs  mb  types. 

Ms  vs  ni[)  values  were  coirected  for  mean  station  magnitude  differences 
prior  to  averaging.  Marshall  and  Basham's  method  was  used  for  MSf 
with  depth  corrections. 


NORTH  LATITUDE 


NORTH  LATITUDE 


37.0i 


35.5—  D  Ms<mb-1.5 

0  mb-1.5<  Ms<  mb-1 


'  mh-.5< Ms 


34.0 


32.5 


31.0 


29.5 


28.0 


26.5 


25.01— 

90.0 


91.5  93.0  94.5  96.0  97.5  99.0  100.5  102. 

EAST  LONGITUDE 

figure  20.  Geographical  distribution  of  various  Ms  vs  me  types. 

vs  mb  values  were  corrected  for  mean  station  magnitude 
differences  prior  to  averaging.  Von  Scggcrn’s  formula  was  used 
for  Ms. 
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figure  21a.  I1  wave  log  A  vs  log  A  plots  before  and  after  the 
application  of  station  corrections  determined  by  the  joint 
magnitude  determination  method. 
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I’igurc  21b,  P  wave  log  A  vs  log  A  plots  before  and  after  the 
application  of  station  corrections  determined  by  the  joint 
magnitude  determination  method. 
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figure  22a.  Rayleigh  wave  log  A  vs  log  A  plots  before  and  after 
the  application  of  station  corrections  determined  by  the  joint 
magnitude  determination  method. 
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figure  22b.  Rayleigh  wave  log  A  vs  log  A  plots  before  and 
after  the  application  of  station  corrections  determined  by 
the  joint  magnitude  determination  method. 


I  i^im  Jo.  Not  of  I'  wave  ami  Kay  I  e 
do  t  o  mi  i  ii  ft  I  by  the  joint  magnitude  d- 


EARTHQUAKE  02/16/1968  5:34:54.2 
33.7  N  95.1  E 
h  =  33  km 

TSIN6HAI  PROVINCE  CHINA 


IR 

i 


EARTHQUAKE  06/28/1968  20:34:55.3 

30.1  N  95.1  E 
h  =  44  km 


Hgurc  24c,  Companion  of  normal  (left)  and  anomalous  event 
seismograms  (right). 
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figure  2U.  Co»pnri*on  of  norw:ii  (left)  and  anonulous  event 
sci$»ograas  (right). 


EARTHQUAKE  06/30/1968  5:05:10.0 
30.2  N  94.8  E 
b  =  42  km 
TIBET 


